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PACS numbers: 12.60.Fr, 12.60.Jv, 13.20.Gd, 95.35.+d, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

Evidence has been accumulated both from astrophysics 
and cosmology that about 1/4 of the energy budget of the 
present universe consists of the so-called (cold) dark mat- 
ter (DM), namely, a component which is non-relativistic 
and neither feels the electromagnetic nor the the strong 
interaction. It is fair to say that the most popular DM 
candidate for a WIMP (Weakly Interacting Massive Par- 
ticle) is the lightest supersymmetric particle (LSP) in su- 
persymmetric models with i?-parity conservation. Leav- 
ing aside the axion and the axino, the superpartners with 
the right properties for playing the role of a WIMP in the 
universe are the gravitino and the lightest neutralino (x) 
- by far, the most discussed case in the literature. 

Although the LHC is running smoothly and collecting 
large amounts of data useful to look for physics beyond 
the Standard Model (SM), other complementary facili- 
ties are certainly needed, especially concerning DM de- 
tection. In fact, DAM A/LIBRA, CoGeNT, and more 
recently CRESST experiments have reported the obser- 
vation of events in excess of the expected background, 
hinting at the existence of a light WIMP [lH3|- How- 
ever, exclusion limits set by other direct searches, such 
as XenonlO [|J and XenonlOO Q, are in tension with the 
above claims. 

In the NMSSM, a light neutralino (as a DM candidate) 
can efficiently annihilate through the resonant s-channel 
via a light pseudoscalar Higgs mediator satisfying the 
requirements from the relic density @. However, fol- 
lowing a scan of the NMSSM parameter space, the au- 
thors of obtained upper limits on the spin- independent 
(SI) x -rLU cleon cross section which are substantially be- 
low the requirements of DAMA and CoGeNT. The spin- 
dependent (SD) cross section (via Z-exchange) was also 
found several orders of magnitude below current exper- 
imental bounds. On the other hand, the authors of @ 
were able to achieve a somewhat larger SI cross section in 
a similar scenario. Admittedly, such cross sections can be 
further enhanced by increasing the s-quark content of the 
nucleon, but the agreement with the low range of DAMA 
results turns out to be only marginally acceptable. 

In this work we revisit SD x-nucleon scattering via 
pseudocalar-exchange in the NMSSM, usually neglected 



in most analyses [7[, which however might be enhanced 
due to a non-perturbative mechanism as later argued. 

II. A LIGHT NMSSM CP-ODD HIGGS BOSON 

The Higgs sector of the NMSSM contains six inde- 
pendent parameters: A, k, A\, A k , tan/3 and ^ e ff, 
whose definitions can be found elsewhere [t| . Notice that 
Me// = As is generated as the vev of the singlet field 
s = (S); it is also useful to define -B e // = Ax + KS - 

Two physical pseudoscalar states appear in the spec- 
trum of the NMSSM as superpositions of the MSSM-like 
state Amssm and the singlet-like state As- In particular 
for the lightest CP-odd Higgs boson 

Ai = cos 9 aAmssm + sm9 A A s , (1) 

where 8a stands for the mixing angle [10]. The A\ 
reduced coupling Xd to down-type quarks and leptons 
(normalized with respect to the coupling of the CP-even 
Higgs boson of the SM) reads 

Xv (A\ — 2ks) 

X c i = cos 8a tan B ~ n x tan 8 , (2) 

M\ + 3kA k s H 1 K ' 

where M\ = 2(j, e f fB e f f / sin 2/3. 

An analysis of a particular region of the NMSSM 
parameter space where Xd can be relatively large at 
high tan/3, together with a light CP-odd Higgs boson 
(rriA 1 ~ 0(10) GeV), was carried out in [l(| although 
without including the relic abundance constraint. Let 
us stress here again that this scenario is quite different 
from the PQ-symmetry-limit (k — > 0) or R-symmetry- 
limit {A K , A\ — > 0), where Xd remains moderate even 
in the large tan/3 limit. Although not motivated by 
any symmetry as in the latter cases, we remark that 
tan/3 ~ l/\[i e ffB e ff \ for large values of tan/3 [llj, giv- 
ing consistency to our scenario which implies (relatively) 
small values of \B e ff\, 

As pointed out in Refs. [12, El a large Xd could induce 
a non- negligible mixing of the CP-odd state and ^(nS) 
hadronic resonances. For the sake of simplicity, we only 
consider the A\ mixing with the nearest (in mass) pseu- 
doscalar resonance, generically denoted hereafter as rjbo- 
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TABLE I: Phase-space corrected leptonic branching fractions, B(T(nS) — > It) (in %), and error bars (summed in quadrature) 
of T(IS'), T(2S), and T(3S) resonances [2Q]. Error bars of the ratios R T /i{nS) are likely overestimated because of expected 
correlations between the numerator and denominator experimental uncertainties. 





B (e+e") 




B(t+t-) 


R T / e (nS) 




T(15) 


2.48 ±0.07 


2.48 ± 0.05 


2.62 ± 0.10 


0.057 ±0.050 


0.057 ±0.046 


T{2S) 


1.91 ±0.16 


1.93 ±0.17 


2.01 ± 0.21 


0.052 ±0.141 


0.041 ± 0.141 


T(3S) 


2.18 ±0.21 


2.18 ±0.21 


2.30 ± 0.30 


0.056 ±0.171 


0.056 ±0.171 


i>(2S) 


0.773 ±0.017 


0.77 ±0.08 


0.772 ±0.100 


0.001 ±0.100 


0.002 ± 0.100 



Thus, the A\ and rjb physical states can be written ap- 
proximately as 



A\ = cos a Aio 
r] b = cos a rjbo 



sin a ijbo 
sin a Aiq 



(3) 
(4) 



where subindex zero refers to unmixed states throughout. 
(The dominant components may of course be reversed if 
a > 7r/4.) In any event, one should keep in mind that 
the Aiq can mix (to a greater or lesser extent) with more 
than a single pseudoscalar hadronic state (see Ref. fisj). 

The strength of the mixing is determined by the angle 
a given by 0, [H| 



sin 2a 



(m 2 



A 



<o) 2 



4 5m 4 



-1/2 



(5) 



where the imaginary part has been neglected and Sm 2 
can be computed by means of a non-relativistic quark 
potential model: 



8m 2 



I 47TO 2 J 



1/2 



\Rr, bo (0)\ X X d , 



(6) 



where v = 246 GeV and R Vbo (0) stands for the radial 
wave function at the origin of the corresponding rjbo state 
(for more details see Ref. 10|). 

As emphasized in Ref. 15], a substantial mixing of a 
O(10) GeV CP-odd Higgs boson with ?/ho resonances can 
modify (hindering) a signal based on direct observation 
of a monochromatic peak in the photon spectrum of ra- 
diative T decays On the other hand, a light CP-odd 
Higgs could still show up as a slight breaking of lepton 
universality in the ratio B TT /Bu ~ 1, where B TT denotes 
the tauonic, and Bu the electronic (£ — e) or muonic 
(I = u) branching ratios of the T resonance, respectively 

In view of the greatly improved accuracy of the re- 
cent measurements of the leptonic BFs (and likely so in 
the forthcoming BaBar analysis of T(3S) decays [13), it 
seems advisable to remove the dependence on the final- 
state lepton mass (rrii) dividing the branching fraction 
(BF) by K{xt) = (1 ± 2x*)(l - 4a^) 1 / 2 , which behaves 
as a (smoothly) decreasing function of xi — m\jM\: 
Bu = Bu/K(xi), with I = e,fi, r. Therefore defining 



B TT - B 



cc 



Brr 
B 



= — 1 ; £ — e,fi 



(7) 



the contribution of a pseudoscalar Higgs to the (inclusive) 
decay rate would imply an enhancement of the tauonic 
mode and therefore small but positive values of TZ T /e [l7| . 

The experimental results obtained from the PDG list- 
ing [20| are shown in Table 1. The good agreement of 
the ipi^S) measurements with the SM expectations, to- 
gether with the systematic disagreement of the T family, 
are consistent with a slight enhancement of the tauonic 
decay mode of T resonances versus the electronic and 
muonic decay modes, due to an extra contribution of a 
light pseudoscalar Higgs boson, which couples to down- 
type quarks (at large tan/3) but of negligible effect for 
up- type quarks in this limit. 

Furthermore, unexpected values for the hyperfine split- 
tings in the bottomonium spectrum: AEhy P (nS) — 
m T ( nS ) - m^fns) {n = 1, 2, 3) [Il,[l3[ can be induced by 
the mixing dlU). In particular, AE hyp (lS) = 69.3 ± 2.8 
MeV obtained by BaBar and CLEO using hindered radia- 
tive T(35) decays [2(| appears to be significantly larger 
than expected from perturbative QCD, estimated to be 
42±13 MeV However, the recent Belle measurement 
of the r]b(lS) mass based on the h(lP) — > rib(lS)j decay 
[H leads to AE hyp (lS) = 59.3 ± 3.1 MeV, in very good 
agreement with lattice NRQCD calculations 2:j . Note 
in passing that the discovery of the r]b(2S) state and its 
mass determination would provide a crucial check when 
compared to the lattice prediction for AEh yp (2S) [23j. 



A. NMSSM scan including WMAP bounds 

In order to assess the impact of the Aiq — r]bo mixing on 
DM phenomenology, we first have to review the present 
bounds from B physics, LEP and DM relic abundance. 
To this aim the latest version of NMSSMTools [H] was 
employed to scan the NMSSM parameter space with mi- 
crOMEGAs turned on in the main code. We focus on a 
narrow mass window for A\ , where current experimental 
constraints still permit large Xd values (2f| under exam- 
ination in this paper. 

The following conditions were required to be satisfied: 

i) 10 < m Al < 10.58 GeV (i.e. below BB threshold). 

ii) Relatively large values of Xd, for large tan fj. 

iii) A light neutralino of mass O(10) GeV, as a WIMP 
candidate satisfying the WMAP bounds. 

The selected ranges of NMSSM parameters correspond 
to a particular but motivated scenario 
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FIG. 1: Xd versus the mass of the CP-odd Higgs boson (m^) 
obtained from our scans on a particular region of the NMSSM 
parameter space for a neutralino mass of a) ~ 3 GeV (black 
circle) b) ~ 3.5 GeV (red diamond); c) ~ 4 GeV (magenta 
asterisk) d) ~ 10 GeV (blue plus); e) ~ 12 GeV (green cross). 



set (i e ft ~ 200 GeV and tan/3 = 45 throughout our 
analysis, while A and k run over the range [0.3,0.5], 
with Ma e [400, 550] GeV. In order to get the highest 
possible Xd values compatible with the bounds, we set 

2 — 3 GeV, dropping to 
- 12 GeV. Note that de- 



A K e [-20, -30] GeV for m x -. 

A K 10 GeV for m x ~ 10 

creasing | A K | implies lowering Xd as can be inferred from 
Eq.@ keeping Ma fixed within the above interval. 

In Fig. 1 we plot Xd versus m.A x for different values of 
the neutralino mass, namely m x = 3, 3.5, 4, 10 and 12 
GeV. For m x close to 5 GeV the values of Xd are very 
small as expected on the grounds of a too efficient annihi- 
lation rate of dark matter due to the resonant condition: 
2m x ~ raA x ■ Notice that values of Xd > 8 are possible 
for m x > 10 GeV, being particularly large (Xd > 10) for 
m x < 4 GeV. It is also worth noting that if condition 
Hi) is removed from the scan, the allowed values of Xd 
become higher, notably about m x = 5 since the afore- 
mentioned resonant condition does not apply anymore. 

Let us also remark that the current limit from B facto- 
ries |H on B[T(3S) -)• 7 + invisible] < (0.7-30) x 10~ 6 , 



for 



1/2 



< 7.8 GeV where s- mv denotes the missing 
invariant-mass squared, actually does not impose any 
bound on the x mass (as > 10 GeV). 



III. NEUTRALINO SCATTERING OFF NUCLEI 

Next we address neutralino elastic scattering off nuclei 
based on our analysis of the NMSSM parameter space. 
As is well-known, the SI scattering cross section is en- 
hanced by a coherent factor proportional to the atomic 
number squared A 2 . For SD interactions, the cross sec- 
tion depends on the total spin of the nucleus and is typ- 
ically a factor A 2 smaller. 

The total WIMP-nucleus cross section has contribu- 
tions from both the SI and SD interactions, though one 



contribution is expected to dominate the other depend- 
ing on the target nucleus (e.g. according to the even/odd 
number of protons and neutrons) and the detection tech- 
nique employed in the experiment. The contributions to 
the SI cross section arise in the interaction Lagrangian of 
the WIMP with quarks and gluons of the nucleon from 
scalar and vector couplings whereas the SD part is at- 
tributed to the axial-vector couplings. Pseudoscalar in- 
teraction is usually neglected because of a strong velocity 
and/or momentum transfer suppression. 

Nevertheless, momentum-dependent interactions have 
been put forward [27j - l29j in order to alleviate the tension 
between the DAMA signal and the null results from other 
experiments. In this work we propose that a significant 
Aiq — f]bo mixing could dramatically modify the SD \~ 
nucleon cross section, in analogy with the well-known 
vector-meson- dominance model (VMD) for electron (or 
real photon) scattering off nuclei [30(, where the virtual 
(or real) photon interacts with nucleons via one of its 
hadronic components. In Fig. 2 we depict two graphs 
illustrating how the mixing of a pseudoscalar (left), or a 
photonic mediator (right), with hadronic resonances can 
modify the effective coupling to the nucleon for WIMP 
and electron scattering, respectively. 

Likewise, a similar mechanism could be envisaged for 
the mixing of scalar resonances (e.g. \bo states) and a 
light enough CP-even Higgs boson. However, it is widely 
accepted that, contrary to the CP-odd Higgs, present 
bounds exclude a O(10) GeV scalar boson with relatively 
large couplings to quarks and leptons (see however [3l|). 
In fact, the lightest scalar Higgs state in our NMSSM 
scan has a mass > 110 GeV, whereby the mixing with 
hadronic states would be negligible. 

The axial-vector and pseudoscalar r/boNN couplings 
defined via 



2M N 



N^ b Nrj b0 - ig Vw NNNj 5 Nr)bo (8) 



lead to a (^-suppression factor at the rate level, where q is 
the momentum transfer of the neutralino to the nucleon. 
An additional q 2 factor stems from the AiqxX vertex, 
yielding altogether a (q 2 /M^^) 2 suppression factor in 
the scattering cross section, where Mjv, x stands for the 
target and neutralino mass, respectively. 

The effective coupling to the nucleon of either the 
mixed A\ or rjb state via its hadronic component reads 



9a/nn = smaxg, )boNN 



%{nn = cosaxg nboNN . (9) 



respectively, where the mixing angle a given by Eq. ([5]). 

In turn, the effective coupling of either the mixed Ai 
or rjb state to the neutralino reads 



eff 
g AiXX 



cos a x gA 10XX 



q eff 



smaxg AloX 



(10) 



respectively, as a result of its Higgs component. 

Let us remark that both A\ and rjb are considered as 
physical states in our model, i.e. mass eigenvalues of the 
full Hamiltonian without transitions among their con- 
stituents apart from interactions with external particles. 
Note, however, that the WIMP scattering off nuclei takes 
place at very low momentum transfer, so that the pseu- 
doscalar mediator in Fig. 2 is quite off-shell. To the extent 
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FIG. 2: Left: Contribution of a pseudoscalar mediator (ei- 
ther a mixed A\ or r/i, state) to the neutralino SD scattering 
cross section. Right: Vector-meson-dominance model in the 
electron scattering off nucleons 1301 ] where vector resonances 
(p, w...) can couple to the exchanged photon. 



that the interaction amplitude is not too sensitive to the 
energy difference between on-shell and off-shell states, 
should our model make sense. 

On the other hand, let us note that the Z-boson cou- 
pling to the neutralino (gzxx) should be of the same or- 
der as the Aiq coupling, but the Aio-coupling to strange 
quarks would be tiny (especially for a dominantly singlet- 
like A w ) as compared to the Z-boson coupling to nucle- 
ons (gzNN) 0- However, the effective coupling ([9]) of 
mixed Ai/rjb states to the latter could be significantly 
enhanced via a non-perturbative effect. Indeed, the ex- 
perimental value B\j] c — > pp] ~ 10~ 3 [20| turns out to 
be unexpectedly large in view of the helicity suppression 
resulting in a perturbative framework. Even after in- 
cluding finite mass effects, a discrepancy of three or four 
orders of magnitude still remains with respect to theoret- 
ical calculations [32| . We shall return later to this point 
of paramount importance in our work, after examining 
the spin-dependent x-nucleon cross section mediated by 
a pseudoscalar Higgs versus a Z-boson. 



A. Pseudoscalar Higgs versus Z-boson exchange 



The x-nucleon scattering cross section with A\ / rjb mix- 
ing (cr™^) can be estimated from the ratio r= 



'SD 



1°. 



SD ' 



where a§ D denotes the Z-exchange SD cross section, ac- 
cording to: 



sin 2 (2a) 



9a 10X x 



9z X x 



9y b oNN 
gZNN 



2 r 



90(GeV) 



1 4 



mA 1 



~] 2 



F 2 [m Al ,m rib ,q 2 } 



(11) 



where F[m Al , ™?7 b , Q 2 } takes into account the interference 
effect due to the sign difference in Eqs.( CTTU|) . At vanish- 
ing q 2 and the range of masses considered in this work, 



|/m 



>n, ■ 



one gets F[m Al ,m nb ,0] ~ \m Al 

Actually a summation over all possible mixed states 
should be understood in (Qj}, albeit with variable weight 
depending on both their mixing strength and the effective 
%o (nS')-nucleon coupling. 



Let us first examine the ratio 

-i 2 



9a 10 xx 



9z xx 



N 2 ,N 2 3 4g 2 cos 2 8 w 



W? 3 



N? 4 } 2 



9l 



cos 2 9 A , (12) 



with Nu denoting the different (bino, higgsino ...) com- 
ponents of the lightest neutralino, and the couplings g\ 
and <?2 satisfy: g 2 /g 2 = ta,n 2 8w, where 9w is the Wein- 
berg angle. Bounds from the invisible decay width of the 



-boson imply that |A 2 S 



N? 4 \ 



< 0.11. Moreover, in 



our scan we always find An close to unity (at large tan/3 
where also A 2 3 » A 2 4 ), and cos 2 (6U) G [1CT 2 - 1CT 1 ]; 
hence we can write 



9A 10 



9z X x 



4A 1 2 1 sin 2 6'vK cos 2 6U 



JV 13 



©(ltr 1 -!) . (13) 



Setting m Al = 10 GeV, sin 2 (2a) ~ 10 _1 (typically ex- 
pected from ([5])) and q 2 ~ (100 MeV) 2 as reference val- 
ues, we are led to 



9v b0 NN 
gZNN 



> c(io 3 



(14) 



so that the Ai/r/b exchange channel would become com- 
parable numerically to the Z-exchange, i.e. r = 0(1). 

Naively one would expect the above ratio to be of or- 
der <x s / yfot em (where ct s and a em denote the strong in- 
teraction and electromagnetic coupling strengths, respec- 
tively) which hardly can yield such a large factor. Yet 
a big enhancement of (TT4"|) could be plausible if a non- 
perturbative mechanism contributes to the ?7bo-nucleon 
coupling, as it likely happens in the t] c (1S) decay into a 
pp pair (though at quite larger q 2 ). 



B. Non-perturbative (instanton-induced) effects 

Indeed, an explanation of the large observed rj c (lS) —> 
pp decay rate seems to require a fundamental modifi- 
cation of the perturbative approach to account for this 
decay mode. Different proposals have been put forward 
in terms of a non-perturbative mechanism: mixing of 
the resonance and gluonium states [33[ , instanton effects 
[34| , intermediate meson loop contribution 35 1 or higher 
Fock components [36$ . Despite many uncertainties, it is 



conceivable that a long-distance contribution also affects 
i]bo(nS) resonances. 

Let us focus hereafter on instanton effects, which play a 
fundamental role in understanding the QCD vacuum and 
many other topics related to hadronic physics, especially 
concerning light hadrons (see e.g. [37[ for a review). Nev- 
ertheless, instanton effects can still be relevant in heavy 
quark systems, e.g. for the non-perturbative gluon con- 
densate in charmonium [38j . 

In particular, the authors of Ref. [3!| studied the in- 
stanton contribution to non-perturbative chiral symme- 
try breaking in proton-proton scattering at high energy. 
The same authors later considered this interaction con- 
tributing to the decay of the rj c resonance into a pp pair 
[3~il [40( | . The idea is that the meson resonance anni- 
hilates into two gluons (perturbative part) that are ab- 
sorbed by instantons, which couple to a baryon pair (non- 
perturbative part). 
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In the following we envisage whether such an 
instanton-induced interaction could still affect the rjb res- 
onance coupling to baryons for a momentum transfer 
from q = m, )i>(is) down to q = 100 MeV. If so, the WIMP 
scattering off nuclei might then bear an unexpected re- 
semblance to p-p elastic scattering (at small — t). 

Instanton effects are usually assumed to depend lin- 
early on the instanton density given by [H], El| : 



dn(p) 
dp 



1 



exp 



2vr 



(15) 



where p denotes the instanton size (such that p < l/q). 
Quite in general, instanton effects are expected to be 
weighted by the Euclidean instanton action of the ex- 
ponential factor in Eq. (|15[) . becoming more relevant at 
smaller q 2 (hence larger a s ). 

Now, in order to obtain an estimate of the 77b (IS) —> pp 
decay rate within this framework, we first rescale the per- 
turbative part of the calculation in [34j to the bottomo- 
nium system according to: 



l^(Q)l 2 

I ^ (0)1 2 



x q£KJ x "' 



2 

2k. 



0.3 , 



(16) 



where the wave function of the spin-singlet state can be 
approximated by the corresponding spin-triplet one [45| . 

Turning now to the non-perturbative part of the cal- 
culation, we will assume (as usual) that the instanton- 
induced coupling to nucleons depends on the number of 
instantons n(l/q) relevant in the process, obtained by 
integration of the instanton density (|15j) over p < l/q. 
Note that the resulting value should be controlled by the 
exponential dependence on p around the inverse mass of 
the heavy resonance. Varying the momentum transfer 



to baryons, from q = m^ c to q 



n(l/q) decreases 



roughly by two orders of magnitude, hence lowering the 
rate (which depends on n(l/q) 2 ) by about four orders of 
magnitude. 

Even though phase space favours the t% — > pp decay 
rate with respect to the i] c by a factor ps 1.3, we finally 
conclude, taking into account both perturbative and non- 
perturbative parts, that the partial width T[rjb — > pp] 
induced by instanton effects should be about four orders 
of magnitude (though with large uncertainties) smaller 
than r[77 c — ¥ pp\. The corresponding BF can be obtained 
making use of the central value for the 776 (IS) full width 
T Vb{ls) = 12.4 MeV, recently found by Belle [22T ]. 

In sum, our order-of-magnitude estimate based on 
instanton-induced interaction reads: 



B[r) b -> pp] w 10" 



10" 



(17) 



Yet the ^(lS) can decay into a pp pair at an observ- 
able rate at B (Super) factories j42(. The experimental 
determination of (|17p thus becomes relevant to uncover 
possible non-perturbative effects (the present upper limit 
being 5 x 10~ 4 20]) associated to the r\b state, permitting 
a reliable comparison with perturbative QCD predictions 
because of the heavier bottom mass. 

On the other hand, as already commented, instanton- 
induced effects should become quite more important at 



low momentum transfer because of the exponential in 
Eq. lfTS"]) . Thus, a sizable non-perturbative effect in r\b 
decays into baryons occurring at q 2 ~ m 2 b should be en- 
hanced (actually less suppressed) at the much lower en- 
ergy scale set by the small momentum transfer in WIMP 
scattering off nuclei, q ~ 100 MeV. However, extrapola- 
tion to such low q value from the r\ c or r\b mass requires 
a further and detailed examination [43j . 

Furthermore, our previous caveat concerning the off- 
shcllncss of the pseudoscalar mediator of Fig. 2 is in or- 
der as indeed q 2 << nrf lb - In particular, its hadronic 
component could display in the nucleon-nucleon vertex a 
different behaviour than an on-shell rjb resonance. Sim- 
ilarly, in a non-covariant language, one can invoke the 
time-energy uncertainty principle implying that the time 
for a virtual bb pair is likely too short for the r\b bound 
state to be formed. Therefore the above calculation of 
the perturbative part cannot not straightforwardly be 
applied. Nonetheless, quantum numbers of the virtual 
bb pair should still correspond to a pseudoscalar state 
thereby permitting two gluons to be emitted, ultimately 
leading to instanton-induced (spin-dependent) effective 
interaction of WIMPs with nuclei. 

Conversely, such a mechanism should not significantly 
affect neutralino annihilation into SM particles via a s- 
channel exchange of a CP-odd Higgs boson, for the en- 
ergy scale would be again of order O(10) GeV, thereby 
avoiding an extra tension with indirect detection limits, 
such as cosmic-ray antiprotons 0, [46| . 



IV. SUMMARY 

In this paper, we have considered a particular sce- 
nario within the NMSSM with both a light neutralino 
and a light CP-odd Higgs boson, the latter sizably mix- 
ing with pseudoscalar rjb resonances. Implicit in our work 
is the idea that non-perturbative effects (e.g. instanton- 
induced interaction) may lead to a non-negligible pseu- 
doscalar contribution to the x-nucleus scattering, thereby 
introducing a momentum-dependent form factor in the 
cross section which might be helpful (see e.g. Ref. |47| ) 
to interpret the results of direct DM search experiments, 
with variable sensitivity along the nuclear recoil energy 
range. 

To conclude we stress that an accurate experimental 
test of lepton universality in T decays, the discovery of 
the ?7b(2S) resonance together with the measurements of 
B[rjb{nS) — > pp] at a (Super) B factory [42| could be 
relevant for a better understanding of DM searches and 
related astrophysical questions. 
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